Energy-landscape network approach to the glass transition 
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We study the energy-landscape network of Lennard- Jones clusters as a model of a glass forming 
system. We hnd the stable basins and the first order saddles connecting them, and identify them 
with the network nodes and links, respectively. We analyze the network properties and model 
the system's evolution. Using the model, we explore the system's response to varying cooling 
rates, and reproduce many of the glass transition properties. We also find that the static network 
structure gives rise to a critical temperature where a percolation transition breaks down the space 
of configurations into disconnected components. Finally, we discuss the possibility of studying the 
system mathematically with a trap-model generalized to networks. 
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I. INTRODUCTION 



In recent years much effort was devoted to the under- 
standing of supercooled liquids and structural glasses, 
and, in particular, the structural arrest taking place at 
the glass transition temperature T g [3, 0] . The numeri- 
cal investigation of the dynamics of supercooled liquids 
and glasses is very hard due to the presence, approach- 
ing T g , of this very slow dynamics Q. An appealing 
approach for understanding this complex dynamics is to 
study the properties of the system's "energy landscape" : 
the dynamics of the system is viewed as the motion of 
the "state point" , described by the 3n-coordinates of all 
particles in the multi-dimensional configuration space, or 
landscape, of the potential energy of the system (n is the 
number of particles). The landscape may be partitioned 
into "basins of attraction" , such that local minimization 
of the potential energy maps any point in a basin to the 
same minimum. In recent years it has been shown that 
the topological details of the basins and the paths con- 
necting them are of great importance in determining the 
properties of glassy systems (e.g., 0,0,0,0,0,0]). 

The representation of the landscape by its basins leads 
to a further simplified view of the energy landscape as a 
network, where the nodes are the basins and the links are 
the saddles connecting them. The energy-landscape net- 
work of a Lennard- Jones (LJ) system has been mapped, 
and some of its properties were extracted H3, EH (for 
energy-landscape networks in proteins and spin systems 
see [lj, [H, H3|). However, the influence of the topology 
of the network on the dynamics of the glass transition 
was never studied. Here, we characterize the networks 
of mono- and bi-disperse LJ systems obtained by min- 
imization of the potential energy, and use a dynamical 
model to study properties such as response to cooling. 
The integration of the landscape picture with network 
theory provides an interpretation of the different critical 
temperatures of the glass transition To (Vogel-Tammann- 
Fulcher temperature [H, 0]) and T g , as well as identifica- 
tion of a new critical temperature T p where a second 
order phase transition separates a phase where a finite 



fraction of the configurations are available, and a phase 
with a vanishing number of accessible states. 

A network model for the glass transition was intro- 
duced in 0, 0]. Here, we take advantage of more sophis- 
ticated network analysis tools such as percolation theory. 
In particular, our approach takes into account the hetero- 
geneity in the number of connections of each basin (i.e., 
its degree fc), which was recently shown to be ubiqui- 
tous in nature and crucial for the understanding of many 
networks' properties [HI]. 



II. THE NETWORK'S STATIC PROPERTIES 

We start with a detailed analysis of the static proper- 
ties of the energy-landscape network. We focus on iso- 
lated small LJ systems of two types: (i) Monodisperse 
LJ system (MLJ) with n = 12,14 particles, (ii) Binary 
80/20 LJ mixture (BLJ) with n = 8 + 2 particles. Our 
network reflects the landscape of potential energy (not 
free energy), or in other words, the entropy is not taken 
into account, since we assume that all basins are equiv- 
alent in terms of the number of internal states they rep- 
resent. We note that the sizes of the systems we study 
are small compared to other systems in which molecular 
dynamics is run [3[. However, this is inevitable since the 
number of nodes increases exponentially with the number 
of particles and thus larger systems are computationally 
much harder to study [101 ]. 

To construct the energy-landscape network, we look for 
basins, the local minima which form the network's nodes, 
and their transition states- first order saddles which con- 
nect two local minima and form the network's links [Hi ]. 
We use the LBFGS algorithm [n| to find the basins, and 
the Eigenmode Method [ijj to find the saddles. Some- 
times more than one first order saddle connect two linked 
basins since the landscape is a high dimensional surface. 
To simplify the network, we consider only the saddle with 
the minimum energy barrier between two linked basins. 
While the BLJ system is known to be glassy [3[ , to avoid 
the crystallization process usually observed in monodis- 
perse systems, we do not consider the state of lowest 



2 



energy when setting up the network [f|. Thus, the two 
systems are expected to be comparable in terms of their 
glassy behavior. 

The MLJ14 network consists of N = 4, 193 nodes 
and M — 58, 628 links. The energies at the nodes are 
distributed approximately normally with mean E(T — > 
00) = —41.5 (Fig. Ha))- The energy barriers 
are distributed approximately exponentially P(AE) = 
where AE = 1.64 is the average energy 
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barrier (Fig. QJb)). As observed in [l0(, we confirm that 
the MLJ14 network is scale-free [l5j . i.e. the degree dis- 
tribution (the probability for a node to have degree k) is 
broad with a tail decaying as P(k) ~ k 1 with 7 « 2.7 
(Fig. LUc)). The energy of a node decreases with its de- 
gree (Fig. [ljd)), meaning that the deepest basins can be 
identified with the network hubs, and are thus accessi- 
ble to/from many other basins [181 ]. The average barrier 
height increases with the degree of the node: AE ~ k e 
with e « 0.43 (Fig . Hl(e )). The average degree of a node's 
nearest neighbors k nn slightly decreases with the node's 
degree (Fig. [ljf )) , meaning hubs have many connections 
to low degree nodes. We also studied "slices" of the net- 
work in the degree (using the g-core method [l9]) and 
energy planes. In both cases, we found that removing 
nodes of low degree, or high potential energy, leaves the 
network connected. 

The MLJ12 system is smaller and contains only N = 
508 nodes and M — 5,407 links, leading to larger fluc- 
tuations in its statistics. Yet the properties of MLJ12 
are qualitatively similar to MLJ14. For MLJ12 we find 
E(T -» 00) = -33.87, AE = 1.16, 7 w 3.1, and e « 0.41 
(Fig. [T]). All the results reported henceforth as ML J are 
for the MLJ14 system, unless explicitly otherwise speci- 
fied. This picture holds true also for the BLJ network, 
with N = 613 nodes_and M = 6, 150 links. Jn_ the BLJ 
network we obtain E(T -> 00) = -30.5, AE = 1.86, 
7 w 3.4, and e « 0.54 (Fig. [I]). 



III. THE NETWORK DYNAMICS 

Next we turn to a characterization of the dynamics of 
the system. We show that application of simple assump- 
tions about the dynamics reproduces many features of 
the glass transition. At high temperatures, kinetic en- 
ergy permits access to most states, while for low tem- 
peratures, mutual access among basins becomes subject 
to considerable activation. In low temperatures near the 
transition (more precisely, below the so-called dynamic 
glass transition temperature Td, where an exponential 
number of meta-stable states appears 0), the dynamics 
is dominated by rare events of collective jumps among 
different stable positions involving many atoms. Thus, 
for low temperatures, we neglect the short time dynam- 
ics which is dominated by small vibrations within the 
basins, and model the dynamics of the system as acti- 
vated jumps between connected states. We assume the 
transition rate between a pair of linked states follows Ar- 
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FIG. 1: Properties of the Lennard- Jones energy-landscape 
network. Shown are results for MLJ14, MLJ12, and BLJ (see 
text). In all figures the data was binned and the average is 
plotted, (a) Distribution of potential energies E of the nodes, 
(b) Distribution of the heights of the energy barriers AE as- 
sociated with the network links, (c) Cumulative distribution 
of node degrees k. Straight lines represent power-law decays 
of the form P(k) ~ k' 1 , with 7 = 2.7,3.1,3.4 for MLJ14, 
MLJ12, and BLJ, respectively, (d) The average potential en- 
ergy of a node E vs. the degree k. (e) The average energy 
barrier to escape from a node AE vs. the node degree k. 
(f) The average degree of node neighbors k nn vs. the node's 
degree k. 



rhenius law: 



Pij 



N -1 



-AEij/T 
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where A-EJy is the height of the barrier separating i and 
j (not necessarily equal to AEji) and the 1/{N — 1) 
factor guarantees that the rate of leaving i, equals to 
Sail links (i j) Pij ls ^ ess t nan 1 f° r an y node. Note there 
is considerable probability for the system to remain at 
the current state. This will turn useful in the character- 
ization of the dynamical slowdown. 

Experiments [l|,l2| and molecular dynamics simulations 
[3j show that supercooling below the melting point results 
in a decrease in the system's energy, up to the temper- 
ature of the glass transition T g . At the transition, the 
system becomes frozen in a disordered configuration, and 
the rate of change of energy with respect to temperature 
decreases abruptly (but continuously) to a value compa- 
rable to that of a crystalline solid. We suggest, that this 
picture, as well as the identification of the glass transi- 



FIG. 2: (Color) MLJio network schematic. The vertical axis 
represents the energy, such that nodes with deeper energy are 
lower (and darker) in the schematic. Links' cusps correspond 
to the energies of the saddles, and nodes' sizes are propor- 
tional to their degree. It can be seen that highly connected 
nodes usually correspond to deeper basins. 



tion temperature T g can be reproduced using our simple 
network dynamics. 

$i (t) , the probability of the system to be at state i at 
time t, evolves according to: 



dt 



1 



N -1 



all links ( 



(2) 

We solve this set of equations numerically by iterating 
Eq. @ once in every time step for the MLJ and BLJ 
networks. We use different cooling rates T(t) = T — Xt, 
where Tj is the initial temperature, and A is the cooling 
rate. We then calculate E(T) = £\ $i(T(t))Ei, where 
Ei is the energy of node i. For infinitely slow cooling, 
the system can be assumed to be in equilibrium, such 
that vanishes for all i. E(T) is calculated by setting 
the Boltzmann distribution $j(T) = e~ Ei l T /Z, where 
Z = J2 t e ~ El/T - The results for BLJ, with A = 0.01 x 
{1/8, 1/16, ...,1/512,0}, are plotted in Fig. Efta). Indeed 
we find that our approach qualitatively reproduces the 
glass-forming behavior. 

We then calculate the heat capacity c = dE/dT (Fig. 
|31Jb)). We associate the temperature for which the heat 
capacity is maximal with the glass transition temperature 
T g . We note that while this association isplausible, it 
cannot be made rigorous. As expected 0, H, H, |2(J, T g 
decreases as the cooling rate becomes slower, approaching 
its equilibrium value T° = 0.67(±0.01) for BLJ (Fig. 
[3](b) , inset). This value of T g is a little higher than the 
glass transition temperature in a large BLJ system (« 
0.45) 3. For MLJ, the picture is similar with T g ° = 
0.47(±0.01). 

Although in our model, microscopic transition rates 
follow Arrhenius law, we show below that the global re- 
laxation times deviate from Arrhenius behavior at low 
temperatures, suggesting that LJ glass forming systems 
are fragile [l[. A global relaxation time is not naturally 
defined for the network. However, we note that as the 
system evolves in time, it explores the phase space in a 
random fashion, according to the transition probabilities 



FIG. 3: (a) Super-cooling in BLJ. The average energy of 
the system E is plotted vs. the temperature T — T — Xt. 
T = 2 and the cooling rates are (top to bottom): A = 
0.01 x {1/8, 1/16, 1/512, 0}, where in each time step we 
iterate Eq. ([2]) once. Zero cooling rate corresponds to the 
equilibrium Boltzmann distribution. At t — we assumed all 
states are equally probable. Similar results are found for MLJ 
(not shown), (b) The heat capacity c = dE/dT. A (bottom to 
top) is same as in (a). Inset: the glass transition temperature 
T g as a function of the cooling rate A. The horizontal line 
corresponds to T° . 
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FIG. 4: Dynamical properties of the energy-landscape net- 
work, (a) The mean first passage time (averaged over all 
sources and destinations) as a function of the inverse tem- 
perature 1/T, for MLJ12 and BLJ. A super- Arrhenius be- 
havior is observed (the slope of the curve increases with 
1/T), suggesting that the system is fragile. The lines are 
fits to Vogel-Tammann-Fulcher law. (b) For several temper- 
atures (top to bottom: T = {2.1, 1.8, 0.3}), the time evo- 
lution of the average energy E(T) was calculated. The y-axis 
shows E(t) — E(t — > 00) (symbols), such that all curves ap- 
proach zero. Curves were fit with a stretched exponential 
E(t) - E(t -> 00) oc exp [~(t/r) p ] (lines). (3 « 0.8 and r 
is between [50,700], and increasing with 1/T. The picture is 
similar for MLJ (not shown). 



given in Eq. (JlJ. Thus, we associate the global relax- 
ation time with the time it takes a random walker with 
transition probabilities as in Eq. JT]) starting at node i, 
to arrive to node j (the first passage time [21|), where 
i and j are randomly chosen, uniformly out of all nodes 
[22| . Given the network and the energy barrier heights, 
the average first passage time can be calculated analyt- 
ically [21]]. In Fig. HJa), we plot the mean first passage 
time (averaged over randomly selected sources and desti- 
nations) as a function of the inverse temperature for BLJ 
and MLJ. A super-Arrhenius behavior is evident, classi- 
fying these systems as a fragile glass [l| . The data seem to 
fit to Vogel-Tammann-Fulcher law r oc exp [A/(T — To)] 
with To ~ 0.1. However, the precise value of To highly 
depends on the simulation details. 

Time dependent quantities can also be studied with 
the network. For example, the evolution of the aver- 
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FIG. 5: A schematic of network percolation, (a) The origi- 
nal network, (b) A fraction q = 1/3 (5/15) of the links are 
removed from the network, (c) The network after removal 
consists of one large cluster of 10 nodes and 2 small clusters 
of one node each. 

age energy of the system at fixed temperature can be 
calculated. We use Eq. ^ and assume that initially 
all states are equally probable. The results are pre- 
sented in Fig. QJb). For short times (up to about 10 3 
time steps) the decay fits to a stretched exponential, 
E(t) -E(t -> oo) = Acxp[-(t/r)' 3 ] with /3 « 0.8 < 1 
0. For longer times (not shown), the decay is exponen- 
tial. As in [f|, the very fast relaxation, corresponding 
to transitions within a basin, is not represented in our 
model. 



IV. PERCOLATION 

In addition to dynamical properties, the network topol- 
ogy gives rise to a static critical temperature T p , where 
the phase space of configurations breaks into discon- 
nected components. This is revealed by percolation the- 
ory applied to the energy-landscape network [23[. Per- 
colation theory is a powerful framework for the study of 
transport in disordered systems. In its simplest form, it 
is engaged in the study of conduction in a lattice in which 
only a fraction p of the sites, or bonds, are conducting 
[24l [25. 26] . This problem is relevant in various contexts 
in which critical phenomena take place, from supercon- 
ductors and gelation to forest fires and oil searching. The 
theory predicts the value of a critical fraction p c above 
which the bulk sample is conducting, as well as the size, 
dimension, total conductance, diffusion coefficients, and 
other properties of the percolation clusters. 

In recent years, percolation theory has been success- 
fully applied to networks to derive criteria for network 
stability. In a percolation process over a network, a frac- 
tion q = 1 — p of the network links is removed [2?], HH . 
A percolation transition occurs when a critical fraction 
q c = 1 — p c of the links is removed such that the net- 
work disintegrates. The critical point where the network 
breaks down is identified by a vanishing size of the largest 
connected cluster as well as a divergence in the size of the 
second largest cluster [24| (Fig- E]). 

At least pictorially, the evolution of the connectivity 



of the energy landscape as the temperature is lowered 
resembles a percolation process. At high enough temper- 
atures, the system has sufficient thermal energy to cross 
most energy barriers. Thus, connected basins are acces- 
sible from each other and the network is intact. At low 
temperatures, links which are associated with a barrier 
of height AE 3> T can practically not be crossed and 
thus can be considered as absent. Thus, as the temper- 
ature is lowered, the network becomes less and less con- 
nected, until reaching the percolation threshold where it 
fully disintegrates. At that point, the system is frozen in 
an isolated region of the landscape, whose size is a zero 
fraction of the entire phase space. A percolation transi- 
tion of the phase space has been predicted long ago for 
spin glasses [231 ] . Here, we use the network representa- 
tion of LJ clusters to show explicitly how the percolation 
transition is realized. 

Since the probability for a link to be "active" decreases 
with decreasing temperature, we suggest that links are 
excluded with probability 1 — e~ AS / T , where AE is the 
link's barrier energy. This way, for high T, all links re- 
main and the network is connected, while for low T many 
links are removed. We then measure (Fig. O the size of 
the largest and second largest cluster (where we define a 
cluster as a set of nodes mutually accessible from each 
other) for ML J and BLJ. A percolation transition is evi- 
dent at T p = 0.26 ±0.01 for MLJ and T p = 0.47 ±0.01 for 
BLJ, indicating a second order phase transition between 
a phase where many configurations are available and a 
phase with a vanishing number of accessible states. 

The percolation transition at T p is expected to take 
place at the final stages of the glass transition, when bar- 
riers become almost impossible to cross, and the system 
freezes in the glassy state. Roughly speaking, the per- 
colation transition temperature T p could be associated 
with the Kauzmann temperature Tk- At Tk, the config- 
urational entropies of the glass and the crystal are equal 
(had the glass transition not intervened), and therefore 
the system is bound to a single, non-crystalline, ideal 
glass state [J Hj]. Similarly, at T p , the system is bound 
to a region of vanishing size of the phase space. In a 
sense, this region in phase space corresponds to the ideal 
glass state in which the system is found at the Kauzmann 
temperature Tk- However we emphasize that this cor- 
respondence is merely descriptive and cannot be made 
more precise. 



V. DISCUSSION 

The understanding of the nature of the glass transi- 
tion is a formidable task, particularly since molecular dy- 
namics cannot approach low enough temperatures. Thus, 
simplified models which capture the essential properties 
of the phenomena are of great value. Representation of 
the multidimensional energy surface as a network is a 
particularly appealing approach, due to recently devel- 
oped network analysis tools. We applied this concept 
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FIG. 6: Percolation transition in Lennard- Jones energy land- 
scape, (a) For MLJ system, we plot the average size of the 
largest and second largest cluster after the removal of each 
link with probability 1 — e _AE//T . Clusters are strongly con- 
nected (i.e., each node can be reached from any other node in 
the cluster). The percolation transition takes place when the 
largest cluster size vanishes and the second largest cluster is 
maximal, (b) Same as (a) for BLJ. 

here, where we studied Lennard- Jones clusters as net- 
works of the stable basins and the links connecting them, 
where each link is associated with an energy barrier. We 
showed that the network approach qualitatively repro- 
duces many properties of the glass transition. It is still 
not known whether quantitative information, such as the 
precise values of T g and other temperatures can also be 
extracted from this kind of analysis. For that purpose, 
larger systems will have to be considered. The similarity, 
in statistical terms, between the networks of n — 12 and 
n = 14 encourages us to believe that similar results, at 
least qualitatively, will be observed in larger systems. 

An alternative approach to circumvent the problem of 
the small system size is a mathematical model which cap- 
tures the main properties of the energy-landscape net- 
work. A naive attempt would be to construct a "general- 
ized trap model" [HI . In a regular trap model the config- 
urations of the system are fully connected, in the sense 
that the system can jump from any state to another. 
However, this is not sufficient to describe the slowing 
down of the dynamics, since it allows transitions which 
do not exist in reality [30j. In a trap model adapted 
to a network, each configuration is linked to precisely k 
other configurations, where it is particularly interesting 
to consider the case of a power-law distribution of de- 



grees which is characteristic of LJ (Section [TTJ) and other 
systems [l2j]. To complete the description of the model, 
one can assume the distribution of energy barriers is ex- 
ponential with mean AE (Section [TTJ) . Despite the at- 
tractiveness of this simple approach, our analysis shows 
[3ll ] that it leads to counterintuitive results. For exam- 
ple, the distribution of time r the system remains in a 
configuration is a power-law P{t) ~ r — (l+W/AJS^ Thus, 
according to model, the typical time the system stays at 
nodes of high degree is small, whereas it expected that 
the system will spend long time at the hubs, since they 
are found at low potential energies (Section [TTJ) . There- 
fore, alternative approaches should be sought for. 

The advantage of the network approach is manifested 
in the application of percolation theory, which provides a 
natural geometrical interpretation of the structural arrest 
taking place at low temperatures [23[ . We studied "bond 
percolation" , where we removed links in which the barrier 
height was high relative to the temperature, to reveal a 
critical temperature where the phase space breaks down 
into small isolated clusters. The study of glassy systems 
with the network approach can be further extended. For 
example, ageing phenomena could be studied, for cither 
the real network or the model, by introducing more com- 
plex correlation functions. Real-space properties such 
as diffusion coefficients and fluctuation-dissipation rela- 
tions could be studied by complementing the network 
with real-space information for each node. In addition, 
similar analysis can be pursued to other systems with 
complex energy landscapes such as proteins (e.g., [321 ]) 
or spin glasses [l4j ]. 
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